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of which has been used to argue that most black-hole growth is obscured' 1 2 1. 
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■^J- . Supermassive black holes underwent periods of exponential growth during 

which we seem them as quasars in the distant Universe. The summed emis- 

sion from these quasars generates the cosmic X-ray background, the spectrum 

Q-f 

There are clear examples of obscured black-hole growth in the form of 'type-2' 

+3 '• 

quasars' 3 but their numbers are fewer than expected from modelling of the 

c3 ■ 

X-ray background. Here we report on the direct detection of a population of 
^ ■ distant type-2 quasars which is at least comparable in size to the well-known 
unobscured type-1 population. We selected objects that have mid-infrared and 
radio emissions characteristic of quasars, but which are faint at near-infrared 
and optical wavelengths. This population is responsible for most of the black 
hole growth in the young Universe and, throughout cosmic history, black-hole 
growth occurs in the dusty, gas-rich centres of active galaxies. 



The population of unobscured type-1 quasars is now well understood out to redshift z ~ 2 
by virtue of optical surveys' 6 ^ which provide direct measurements around the "break" in the 
luminosity function where the quasar population outputs most of its luminosity density. 
Current understanding of the obscured type-2 quasar population is far less complete. 

Examples of type-2 quasars have, for many years, been confined mainly to radio-selected 
samples. Radio emission is unaffected by dust and, at low frequencies, is not strongly de- 
pendent on viewing angle so the simplest interpretation of optical spectroscopic follow-up of 
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radio-selected samples is that luminous quasars divide roughly equally into obscured (type- 
2) and unobscured (type-1) objects' 7 '. However, 'radio-loud' quasars are rare and may be 
atypical: their large-scale radio jets may well modify their environments in ways which in- 
fluence whether or not the nucleus is obscured, and hence the ratio of quasars in the type-1 
and type-2 classes' 8 ! . 

More recently, deep X-ray surveys and optical surveys have also been successful in identi- 
fying type-2 quasars' 5 ' 9 ' 10 '. However, optical spectroscopic follow-up has yet to find a type-2 
population which can account fully for the X-ray background' 1 ' 9 ' 11 '. This is believed to be 
because a substantial fraction of type-2 quasar nuclei are hidden by 'Compton thick' mate- 
rial (with neutral gas column densities N H ^ 10 28 uT 2 ) even when observed at high X-ray 
energies' 1 ' 12 '. 

type-2 quasars (with visual extinction Ay ^ 5 magnitudes towards their nuclei) fail to 
outshine their host galaxy at ultra-violet and optical wavelengths. However, at longer wave- 
lengths the obscuration becomes small so the mid-infrared to radio properties of type-2 
quasars should be similar to those of their type-1 counterparts, unless the obscuration or 
redshift becomes extreme. Indeed a survey at 60 /xm with the Infrared Astronomy Satellite 
(IRAS) led to the discovery of the first high-redshift radio-quiet type-2 quasar: the z = 2.3 
galaxy IRAS F10214+4724' 4 '. However, this object was discovered only because its flux was 
magnified by a factor ~ 50 by a gravitational lens' 13 ' . From the presence of high-excitation 
narrow lines and scattered broad lines' 14 ' it clearly has an obscured type-2 nucleus but its 
X-ray emission seems to be obscured by Compton-thick material' 12 '. 

With the advent of the Spitzer Space Telescope' 15 ' , with a sensitivity > 100 times greater 
than IRAS and > 10 than the Infrared Space Observatory (ISO), type-2 quasars similar 
to IRAS F10214+4724 can be found from their mid-infrared emission without the aid of 
gravitational lenses. We chose the following selection criteria to find high-redshift type-2 
quasars in the Spitzer First Look Survey (FLS): mid-infrared 24-yum flux density S , 24 (Um > 300 
/xJy; near-infrared 3.6-/im flux density S 3 , e ^m ^ 45 /iJy; and 350 /iJy ^ Si.4GHz ^ 2 mJy, 
where Si.4gh z is the 1.4-GHz radio flux density. 

The mid-infrared (5 , 2 4 Mm ) criterion was chosen to obtain a reliable (7-er) catalogue from 
the Spitzer FLS data from the MIPS instrument (Fadda, D., et al, in preparation). At z = 2 
this means the quasar luminosity will be > 0.2 L* uasar (where L* quasar is the break luminosity 
which corresponds to Mb = —25.7 magnitudes' 16 '); this means we can detect sources around 
the break in the luminosity function (see Supplementary Information). Observed 24 /im 
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corresponds to emitted 8 /im at z — 2, and dust extinction is negligible at this wavelength 
unless the obscuring column is extreme. (We adopt a ACDM cosmology with the following 
parameters: h = H /(100 km s^ 1 Mpc' 1 ) = 0.7; Q m = 0.3; Q A = 0.7.) 

The upper limit on the 3.6 /im flux density (obtained using the IRAC instrument l 17 ') 
was chosen to select only the most distant type-2 quasars. At z ^ 2 the detected 3.6-/xm flux 
density corresponds to light emitted at A ^ 1.2 /im, so we are able to eliminate unobscured, 
type-1 quasars whose flux would exceed this limit (FIG 1). Indeed, for type-2 quasars the 
<S3.6|um emission is likely to be dominated by starlight from the host galaxy, allowing us to 
estimate 'photometric redshifts' z phot (see Supplementary Information). Our S 3 ^ m criterion 
corresponds to z phot > 1.4 (see caption to Table 1 for more details). 

The Sssum criterion includes objects without IRAC detections which we expect to have 
the highest redshifts. Since the 24-/im positions are accurate to ~ 1 arcsec, the FLS radio 
positions' 18 ', accurate to ~ 0.5 arcsec, were better for spectrosocopic follow-up. However, 
the main importance of the radio selection criteria is to ensure that the candidates are 
quasars rather than starburst galaxies. We chose a lower limit on Si.4GHz well above the 
level reached by high-redshift (submillimetre-selected) starburst galaxies without the benefit 
of gravitational lensing' 19 ' and an upper limit to filter out the radio-loud objects, whose 
extended jets might complicate interpretation. 

FIG. 1 shows the location of the 21 type-2 quasar candidates that met our selection 
criteria in the versus Sssum plane and Table 1 summarises their properties. The 

candidates were found from a sky area of 3.8 deg 2 in the Spitzer FLS. No type-1 quasar 
(A v = 0) that is bright enough to be detected at 24 jum will be faint enough at 3.6 /im to 
be selected (FIG. 1). "Blind" low- resolution optical spectroscopy ' 20 1 was performed on our 
21 candidates, 10 of which yielded clear type-2 spectra with narrow emission lines giving 
redshifts in the range 1.4<z<4.2. These objects are clearly type-2 quasars either because 
they have high-excitation lines (e.g. inset to FIG. 1) or because the rest-frame equivalent 
widths of their Ly-a lines are > 100 nm and hence significantly larger than those seen 
in starbursts. Of the 11 objects that yielded no redshifts, all but one (which shows faint 
red continuum) have completely blank spectra, showing there is probably no contamination 
from lower-redshift (z< 1.4) starbursts since they would have shown [Oil] line emission, or 
if highly obscured, at least some continuum. We believe these blank-spectrum objects are 
also type-2 quasars with z>1.4 and there is no compelling argument against them having 
z > 2 since the Lyman-a line can easily be extinguished if the host galaxy is dusty on large 
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(kpc) scales. We are probably seeing two types of type-2: the objects with narrow emission 
lines are obscured by the torus' 3 ', while the blank-spectrum objects are probably obscured 
by a starbursting host galaxy and we are seeing the coeval growth of supermassive black 
hole and host galaxy. Growing supermassive black holes embedded in a starburst have been 
found at z ~ 2 from X-ray measurements of submillimetre-selected galaxies' 21 '. 

To interpret our results in terms of the 'quasar fraction' q - the ratio of the number 
of type-1 quasars to the total number of type-1 and type-2 quasars - we need to predict 
the average number < JVi > of type-1 quasars meeting identical S^/an and Si.4ghz selection 
criteria, and having matched redshift and sky area selection functions. We estimated a 
probability distribution p(< Ni >) to account for the various uncertainties, the gaussian 
approximation of which (4.3l 2 ; 2 ) represents the expected number of type-1 quasars expected 
to adhere to both our 24-//m and 1.4-GHz selection criteria' 22 ', and with z ^ 2 in a 3.8 deg 2 
patch; this number would be 15-times higher without any radio selection criteria. 

FIG. 2 shows the resulting conditional probability distribution of the quasar fraction q 
given our new data on type-2 quasars and our background knowledge of the type-1 quasar 
population, namely p(g|data, {type — 1 quasar}). This bayesian approach allows us to ac- 
count properly for small number statistics (see FIG. 2 and Supplementary Information). It 
is difficult to judge whether the blank-spectrum objects with photometric redshifts z p hot > 2 
should be included (dashed red line in FIG. 2) until we know more about the dust content of 
their host galaxies, but considering just the objects with spectroscopic redshifts (solid blue 
line in FIG. 2) is clearly a conservative estimate. 

The most likely systematic problem with our analysis is that we have implicitly assumed 
identical radio properties' 22 ' for type-1 and type-2 quasars matched in intrinsic luminos- 
ity. However, 'unified models'' 3 ' predict Doppler-boosted radio emission from the weak jets 
aligned closer to the line-of-sight for type-1 quasars' 23 ' and hence, on average, lower radio 
luminosities for type-2 quasars. This would yield lower derived values of q and strengthen 
our conclusion that most quasars are obscured. It is of course possible to postulate that the 
optical-to-radio correlation evolves with redshift, but this is currently an unknown (see also 
Supplementary Information). 

The median redshift of our spectroscopic sample is at z = 2 which is exactly where the 
optical luminosity density of the type-1 quasar population is at its peak' 6 '. Our sample shows 
that during this 'epoch of quasar activity' we can be confident that about half of all high- 
luminosity quasars are obscured and if some of the blank-spectrum objects turn out to have 
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z > 2 then it is likely that type-2 quasars outnumber type-1 quasars by a significant factor. 
It is also possible that we are missing type-2 quasars with an extreme obscuring column 
(Ay > 100) which could be due, for example, to the obscuring torus being completely edge- 
on to the line of sight, but this would only strengthen our result, which has been predicted 
both by models for the X-ray background' 1 ', and by models predicting the surface density of 
faint AGN in relatively small sky area surveys' 24 '. This is, however, the first direct detection 
of a high-redshift type-2 quasar population likely to outnumber the type-1 population. We 
note that the solid blue line (which compares narrow-line or 'torus-obscured' type-2s with 
type-Is) is consistent with the quasar fraction from radio-loud samples' 7 ', while the dashed 
red line (which includes both 'torus-' and 'host-obscured' type-2s) has enough obscured 
objects to account for the X-ray background' 1 ' . 

FIG. 2 shows the expected ratio for two 'receding torus'' 25 ' 26 ' 27 ' models in which the 
solid angle subtended by the obscuring torus decreases (so that q increases) systematically 
with increasing quasar luminosity. Since we have suggested that some of the blank-spectrum 
objects are obscured by kpc-scale dust, it is only the ratio of spectroscopically-confirmed 
type-2 quasars to type-1 quasars (the solid blue line) that is relevant to the model predictions, 
thus these models are not excluded. 

To conclude, we consider our results alongside those of the X-ray surveys at lower 
redshifts' 9 ' (z < 2) which find that most accretion at low redshift occurs in type-2 quasars 
and type-2 Seyferts (their lower luminosity analogues). We deduce that throughout cosmic 
history, black hole growth seems to have been concentrated in obscured regions. This is in 
good agreement with predictions from the X-ray background' 1 ' and implies, from compar- 
isons between the integrated luminosity density of quasars (both type-1 and type-2) and the 
local space density of relic black holes' 28 ', that black hole growth occurs in short, efficient 
spurts in the cores of forming galaxies' 2 ' . 
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Table 1. Basic data on the 21 type-2 quasars in our sample. The J2000.0 positions are from the Spitzer FLS radio catalogue' 18 ' . 
RA, right ascension, Dec, declination. The MIPS 24-fj.m flux density S24/jm is obtained by point spread function (PSF) fitting! 29 ' 
as all objects are point sources at the ~ 6 arcsccond resolution of the MIPS observations, with positional errors ~ 1 arcsecond. 
It has a typical error of ±10 — 15 per cent. The 1.4-GHz flux density is the peak value (in /ijy per beam) from the radio 
catalogue! 18 ' . The IRAC 3.6- and 4. 5-^tm flux densities are measured in 5-arcsec-diameter apertures! 17 ' and have typical errors 
of ±10 per cent. The photometric redshifts arc calculated by assuming that the 3.6-^tm flux density is dominated by starlight 
from a 2L* aJ galaxy (L* al is the break in the galaxy luminosity function, following the models of ref. ! 30 ' and passive evolution). 
The criterion S3Sfim ^ 45 /ijy corresponds to a limiting photometric redshift z p hot ^ 1-4; this was chosen to allow for scatter 
in the photometric redshift estimation, whilst still filtering out type-1 quasars and low-redshift contaminants like radio galaxies. 
All of the 21 candidates were observed spectroscopically with the dual-beam ISIS instrument at the William Hcrschcl Telescope, 
on July 2004 and April 2005. 'Blind' low-resolution spectroscopy' 20 ' was performed by offsetting from nearby stars to the radio 
positions, with exposure times of ~ 30 minutes and a continuous wavelength coverage across the entire visible band. Objects 
with a 1 have only one definite line in their optical spectra, which is taken to be Ly-a by virtue of extreme equivalent width, 
blue-absorbed line profile and lack of other lines supportive of alternative identifications. 
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Figure 1. Diagrammatic representation of the infrared selection criteria. Logio(24 fim flux density 524pm / A*Jy) plotted 
against logio(3.6 fim flux density £3.6^771/ A*Jy)- The dotted lines represent the flux density cuts of our selection criteria. Note 
that all objects plotted needed to meet the criteria on radio flux density described in the text, and that the candidates were 
originally selected from the preliminary 3.6 and 24 fim catalogues, but have had their fluxes slightly revised. This explains 
why two candidates are slightly outside the selection boundaries (see Supplementary Information). The parent population is 
plotted as orange crosses, or arrows when 524pm is an upper limit. Our 21 candidates are plotted in the top left corner, with 
black triangles for candidates with a spectroscopic redshift and black diamonds for candidates without one. Horizontal arrows 
represent upper limits for 53.6pm- None of the candidates showed type-1 quasar spectra. The spectrum of the z = 4.174 object 
is shown in the inset. The magenta crosses show 524pm and 53.6pm for a 2L* ai elliptical galaxy at z = 1 and z = 2; starlight is 
assumed to evolve passively following a stellar population synthesis model t 30 l . The other colours represent an L* uasar quasar 
(spectral energy distribution from refl 31 !) with varying amounts of visual extinction Ay applied to 524pm and S3.6p m and 
plotted on the curves (assuming a Milky- Way-like extinction curve from refJ 32 '); these are plotted at two different redshifts, 
and the quasar population is assumed to undergo pure luminosity evolution' 16 ! . At all redshifts, any type-1 (Ay = 0) quasar 
bright enough to make the 524pm > 300 /ijy criterion will also be too bright at 3.6-fim to make it into our sample. Inset, the 
optical spectrum of AMS16, showing Lyman-c?, N V and C IV emission lines. 
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Figure 2. Probability distributions for the quasar fraction. The normalised posterior probability p(g|data, {type — 1 quasar}) 
of the quasar fraction q given our new data on type-2 quasars, and background information on type-1 quasars. The two 
lines represent different interpretations of the number of type-2 quasars at 2 ^ 2 determined from our new data. To ob- 
tain p(g|data, {type — 1 quasar}) we used Bayes' theorem with a uniform prior over the range q 1 so that the 
p(g|data, {type — 1 quasar}) = p(data|g, {type — 1 quasar}). We then calculated this likelihood function, at each q, as a 
Poisson distribution for the observed (integer) number of type-2 quasars at z 2 where the mean number expected 
< N2 >= (1 — q) < Ni > /q. Because of uncertainties in the background information on < Nj >, we had to evaluate 
the likelihood at each q and < N\ > and then integrate p< N\ > X the likelihood over < Ni >. The solid blue line is the 
posterior probability if we believe the only candidates with z 2 are the 5 whose spectroscopic redshifts confirm this. If, 
additionally, we use the photometric redshifts of the candidates with blank spectra (6 additional candidates with z ^ 2), then 
we obtain the red curve. The best estimate for q is probably somewhere between the two curves. The two vertical lines represent 
the predictions of receding torus models, one with a fixed torus height (marked 'RTl')l 26 ' and one ('RT2') with a torus height 
varying with luminosity! 27 ) . 



